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Abstract. Dark matter has been introduced to explain substantial mass deficits 
noted at different astronomical scales, in galaxies, groups of galaxies, clusters, 
superclusters and even across the full horizon. Dark matter does not interact with 
baryonic matter except gravitationally, and therefore its effects are sensed only on 
the largest scales. Although it is still unknown whether dark matter consists of 
particles or of a field or has some other nature, it has a rich phenomenology. This 
review summarizes all the astrophysical and cosmological probes that have produced 
overwhelming evidence for its existence. The breadth of the subject does not permit 
details on the observational methods (the reference list then helps), thus the review is 
intended to be useful mainly to cosmologists searching to model dark matter. 
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1. Introduction 

Apparently the matter content of the Universe is dominated by an unknown form of 
dark matter (DM) without interactions with ordinary baryonic matter, perhaps not 
even with itself. It only interacts via the gravitational field, manifesting its effects on 
astrophysical and cosmological scales. The purpose of this review is to summarize the 
phenomenology of all such effects, that can serve as probes of dark matter. Regardless 
of the ultimate, correct explanation of its particle nature or field nature, theory needs 
to address all these effects. 

This review does not cover the historical development, except by glimpses, because 
the rapid development of observational means tends to render all discoveries older than 
a decade unimportant. 

Beginning from the first controversial conclusions from the motion of stars near the 
Galactic disk on missing matter in the Galactic disk (Sec. 2), and that of Fritz Zwicky 
in 1933 [1] of missing matter in the Coma cluster (Sec. 3), we describe the kinematics 
of virially bound systems (Sec. 3) and rotating spiral galaxies (Sec. 4). An increasingly 
important method to determine the weights of galaxies, clusters and gravitational fields 
at large, independently of electromagnetic radiation, is lensing, strong as well as weak 
(Sec. 5). Next follows a discussion of dark matter in elliptical galaxies (Sec. 6) and 
mass-to-light ratios which probe dark matter in all systems, notably in dwarf spheroidals 
(Sec. 7). Different ways to measure missing mass in groups and clusters derive from the 
comparison of visible light and X-rays (Sec. 8). Mass autocorrelation functions relate 
galaxy masses to dark halo masses (Sec. 9). 
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In radiation the most important tools are the temperature and polarization 
anisotropies in the Cosmic Microwave Background (CMB) (Sec. 10), which give 
information on the mean density of both dark and baryonic matter as well as on 
the geometry of the Universe. The large scale structures of matter exhibit similar 
fluctuations evident in the Baryonic Acoustic Oscillations (BAO) (Sec. 11). The 
amplitude of the temperature variations in the CMB prove, that galaxies could not have 
formed in a purely baryonic Universe (Sec. 12). Simulations of large scale structures 
also show that DM must be present (Sec. 13). The best quantitative estimates of the 
density of DM come from overall parametric fits to cosmological models, notably the 
Cold Dark Matter model 'ACDM' with a cosmological constant A, of CMB data, BAO 
data, and redshifts of supernovae of type la (SNela) (Sec. 14). A particularly impressive 
testimony comes from merging clusters (Sec. 15). We conclude this review with a brief 
summary (Sec. 16). 

2. Stars near the Galactic disk 

In 1922 the Dutch astronomer Jacobus Kapteyn [2] studied the vertical motions of all 
known stars near the Galactic plane and used these data to calculate the acceleration 
of matter. This amounts to treating the stars as members of a "star atmosphere", a 
statistical ensemble in which the density of stars and their velocity dispersion defines a 
"temperature" from which one obtains the gravitational potential. This is analogous to 
how one obtains the gravitational potential of the Earth from a study of the atmosphere. 
Kapteyn found that the spatial density is sufficient to explain the vertical motions. 

Later in the same year the British astronomer James Jeans [3] reanalyzed Kapteyn's 
data and found a mass deficit: to each bright star two dark stars had to be present. 
The result contradicted grossly the expectations: if the potential provided by the known 
stars was not sufficient to keep the stars bound to the Galactic disk, the Galaxy should 
rapidly be losing stars. Since the Galaxy appeared to be stable there had to be some 
missing matter near the Galactic plane. 

In 1932 the Dutch astronomer Jan Hendrik Oort [3] reanalyzed the vertical motions 
and came to the same conclusion as Jeans. There was indeed a mass deficit which Oort 
proposed to indicate the presence of some dark matter in our Galaxy. The possibility 
that this missing matter would be nonbaryonic could not even be thought of at that time. 
Note that the first neutral baryon, the neutron, was discovered by James Chadwick [5] 
only in the same year, in 1932. 

However, it is nowadays considered, that this does not prove the existence of DM 
in the disk. The potential in which the stars are moving is not only due to the disk, but 
rather to the totality of matter in the Galaxy which is dominated by the Galactic halo. 
The advent of much more precise data in 1998 led Holmberg & Flynn [6] to conclude 
that no DM was present in the disk. 

Oort determined the mass of the Galaxy to be 10 11 M sun , and thought that the 
nonluminous component was mainly gas. Still in 1969 he thought that intergalactic gas 
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made up a large fraction of the mass of the universe [7] . The general recognition of the 
missing matter as a possibly new type of non-baryonic DM dates to the early eighties. 

3. Virially bound systems 

The planets move around the Sun along their orbits with orbital velocities balanced by 
the total gravity of the Solar system. Similarly, stars move in galaxies in orbits with 
orbital velocities v determined by the gravitational field of the galaxy, or they move 
with velocity dispersion a. Galaxies in turn move with velocity dispersion a under the 
influence of the gravitational field of their environment, which may be a galaxy group, 
a cluster or a supercluster. In the simplest dynamical framework one treats massive 
systems (galaxies, groups and clusters) as statistically steady, spherical, self-gravitating 
systems of N objects with average mass m and average velocity v or velocity dispersion 
a. The total kinetic energy E of such a system is then (we now use a rather than v) 

E = (1/2)W. (1) 

If the average separation is r, the potential energy of iV(iV — l)/2 pairings is 

U = -(1/2)N(N- l)Gm 2 /r . (2) 
The virial theorem states that for such a system 

E = -U/2 . (3) 
The total dynamic mass M^ yn can then be estimated from a and r 

M dyn = Nm = 2ra 2 /G . (4) 
This can also be written 

a 2 oc {M dyn /L)IR , (5) 

where I is a surface luminosity, R is a scale, and Md yn /L is the mass-to-light ratio. 
Choosing the scale to be the half light radius R e , this implies a relationship between the 
observed central velocity dispersion ctq, I e and R e called the Fundamental Plane, of the 
form 

R e oc K) a (/ e ) fe • (6) 

The virial theorem predicts the values a = 2, b = 1 for the coefficients. This relationship 
is found in ellipticals E] and in some other types of stellar populations, but with 
somewhat different coefficients. 

3. 1 . Halo density profiles 

The shapes of DM halos in galaxies and clusters need to be simulated or fitted by 
empirical formulae. Mostly the shape is taken to be spherically symmetric so that the 
total gravitating mass profile M(r) depends on three parameters: the mass proportion 
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in stars, the halo mass and the length scale. A frequently used radial density profile 
parametrization is 

Pdm{t) = Po /[(r/r s r (1 + r/r s ) 3 " a ] , (7) 

where po is a normalization constant and < a < 3/2. Standard choices are a = 1 for 
the Navarro- Frenk- White profile (NFW) [10J, and a = 3/2 for the profile of Moore & 
al. [IT], both cusped at r = 0. 

Another parametrization is the Einasto profile ([12] and earlier references therein) 

Pdm{t) = p e ex V {-d n [(r/r e )^ n - 1]}, (8) 

where the term d n is a function of n such that p e is the density at r e , which defines a 
volume containing half of the total mass. At r = the density is then finite and cored. 
The Burkert profile [13J has a constant density core 

p DM (r) = p /[(l + r/r s )(l + (r/r s ) 2 )] , (9) 

which fitted dwarf galaxy halos well in 1995, but no longer does so, see Sec. 7. 

Some clusters are not well fitted by any spherical approximation. The halo may 
exhibit a strong ellipticity or triaxiality in which case none of the above profiles is good. 

The dependence of the physical size of clusters on the mass, characterized by the 
mass concentration index c = r V i r /r s , has been studied in ACDM simulations [14]. At 
intermediate radii c is a crucial quantity in determining the density shape. 




Figure 1. Density profile of matter components enclosed within a given radius r in 
the Coma cluster, versus r/r V i r . From E. L. Lokas & G. A. Mamon 15j . 



Astrophysical and cosmological probes of dark matter 



6 



3.2. The Coma cluster 

Historically, the first observation of dark matter in an object at a cosmological distance 
was made by Fritz Zwicky in 1933 pQ. While measuring radial velocity dispersions of 
member galaxies in the Coma cluster (that contains some 1000 galaxies), and the cluster 
radius from the volume they occupy, Zwicky was the first to use the virial theorem to 
infer the existence of unseen matter. He found to his surprise that the dispersions 
were almost a factor of ten larger than expected from the summed mass of all visually 
observed galaxies in the Coma. He concluded that in order to hold galaxies together the 
cluster must contain huge amounts of some non-luminous matter. From the dispersions 
he concluded that the average mass of galaxies within the cluster was about 160 times 
greater than expected from their luminosity (a value revised today), and he proposed 
that most of the missing matter was dark. 

Zwicky's suggestion was not taken seriously at first by the astronomical community 
which Zwicky felt as hostile and prejudicial. Clearly, there was no candidate for the 
dark matter because gas radiating X-rays and dust radiating in the infrared could not 
yet be observed, and non-baryonic matter was unthinkable. Only some forty years later 
when studies of motions of stars within galaxies also implied the presence of a large 
halo of unseen matter extending beyond the visible stars, dark matter became a serious 
possibility. 

Since that time, modern observations have revised our understanding of the 
composition of clusters. Luminous stars represent a very small fraction of a cluster 
mass; in addition there is a baryonic, hot intracluster medium (ICM) visible in the X- 
ray spectrum. Rich clusters typically have more mass in hot gas than in stars; in the 
largest virial systems like the Coma the composition is about 85% DM, 14% ICM, and 
only 1% stars [T5] . 

In modern applications of the virial theorem one also needs to model and 
parametrize the radial distributions of the ICM and the dark matter densities. In 
the outskirts of galaxy clusters the virial radius roughly separates bound galaxies from 
galaxies which may either be infalling or unbound. The virial radius r„j r is conventionally 
defined as the radius within which the mean density is 200 times the background density. 

Matter accretion is in general quite well described within the approximation of the 
Spherical Collapse Model. According to this model, the velocity of the infall motion and 
the matter overdensity are related. Mass profile estimation is thus possible once the 
infall pattern of galaxies is known [16J. 

In Fig. 1 the Coma profile is fitted [15] with Eq. ^ with a = which describes 
a centrally finite profile which is almost flat. The separation of different components in 
the core is not well done with Eq. ([7]) because the Coma has a binary center like many 
other clusters [T7j . 
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Figure 2. Density profile of matter components in the cluster AC 114, enclosed within 
a given projected radius. From M. Sereno & al. (18) . 

3.3. The AC 114 cluster 

Dark matter is usually dissected from baryons in lensing analyses by first fitting the 
lensing features to obtain a map of the total matter distribution and then subtracting 
the gas mass fraction as inferred from X-ray observations [T9l 120] . The total mass map 
can then be obtained with parametric models in which the contribution from cluster- 
sized DM halos is considered together with the main galactic DM halos [21] • Mass in 
stars and in stellar remnants is estimated converting galaxy luminosity to mass assuming 
suitable stellar mass to light ratios. 

One may go one step further by exploiting a parametric model which has three 
kinds of components: cluster-sized DM halos, galaxy-sized (dark plus stellar) matter 
halos, and a cluster-sized gas distribution [T71 [18]. As an example we show the results 
of such an analysis of the dynamically active cluster AC 114 in Fig. 2. 

In systems of merging clusters DM may become spatially segregated from baryonic 
matter and thus observable. We shall meet several such cases in Sec. 15. 

3.4- The Local Group 

The Local Group is a very small virial system, dominated by two large galaxies, the 
M31 or Andromeda galaxy, and the Milky Way. The M31 exhibits blueshift, falling in 
towards us. Evidently our Galaxy and M31 form a bound system together with all or 
most of the minor galaxies in the Local Group. The Local Group extends to about 3 
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Mpc and the velocity dispersions of its members is about 200 km s _1 . 

In this group the two large galaxies dominate the dynamics, so that it is not 
meaningful to define a statistically average pairwise separation between galaxies, nor 
an average mass nor an average orbital velocity. The total kinetic energy E is still given 
by the sum of all the group members, and the potential energy U by the sum of all the 
galaxy pairs, but here the pair formed by the M31 and the Milky Way dominates, and 
the pairings of the smaller members with each other are negligible. 

An interesting recent claim is, that the mass estimate of the Local Group is also 
affected by the accelerated expansion, the "dark energy". A. D. Chernin & al. [22] have 
shown that the potential energy U is reduced in the force field of dark energy, so that 
the virial theorem for N masses with baryocentric radius vectors r-j takes the form 



is a correction which reduces the potential energy due to the background dark energy 
density p v . In the Local Group this correction to the mass appears to be quite 
substantial, of the order of 30% — 50%. 

The dynamical mass of the local group is 3.2 — 3.7 x 10 12 solar masses whereas the 
total visible mass of the Galaxy + M31 is only 2 x 10 11 solar masses. Thus there is a 
large amount of dark matter missing. 

3.5. The local Universe 

In a large volume beyond the local group, Tully in 1984 [23] measured the velocities 
of 2367 galaxies with radial velocities below 3000 km s _1 . He found that the mass 
density parameter (which is normalized to the critical mass) in this "Local Universe" 
was Q m = 0.08, in clear conflict with the global value, fl m , g i bai — 0.27 ± 0.02 (as we 
shall see in Sec. 14). 

More recently Karachentsev [24] has extended this analysis out to a volume of a 
diameter of 96 Mpc, containing 11 000 galaxies appearing single, in pairs, in triplets and 
in groups. Most of them belong to the Local Supercluster and constitute < 15% of the 
mass of Virgo. The radial velocities are v < 3500 km s _1 . These galaxies can be treated 
as a virial system with average density fl m ,iocai = 0.08 ± 0.02, again surprisingly small 
compared to the global density. Karachentsev quotes three proposed explanations for 
this mass deficit. 

- Dark matter in the systems of galaxies extends far beyond their virial radius, so 
that the total mass of a group or cluster is 3 - 4 times larger than the virial estimate. 
However, this contradicts other existing data. 

- The diameter of the considered region of the Local universe, 90 Mpc, does not 
correspond to the true scale of the "homogeneity cell" ; our Galaxy may be located inside 
a giant void sized about 100 - 500 Mpc, where the mean density of matter is 3 to 4 



E = -(1/2)U + U 2 , 



(10) 



where U is defined as in Eq. @, and 

U 2 = -(4irp v /3) Em^r 2 



(11) 
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times lower than the global value. However, the location of our Galaxy is characterized 
by an excess, rather than by a deficiency of local density at all scales up to 45 Mpc. 

- Most of the dark matter in the Universe, or about two thirds of it, is not associated 
with groups and clusters of galaxies, but distributed in the space between them in the 
form of massive dark clumps or as a smooth "ocean". It is as yet difficult to evaluate 
this proposal. 

Clearly the physics in the Local Universe does not prove the existence of dark 
matter, rather it brings in new problems. 

4. Rotation curves of spiral galaxies 

Spiral galaxies are stable gravitationally bound systems in which visible matter is 
composed of stars and interstellar gas. Most of the observable matter is in a relatively 
thin disc, where stars and gas rotate around the galactic center on nearly circular orbits. 
The galaxy kinematics is measured by the Doppler shift of well-known emission lines of 
particular tracers of the gravitational potential: HI, CO and H Q . 

If the circular velocity at radius r is v in a rotating galaxy with mass M(r) inside 
r, the condition for stability is that the centrifugal acceleration v/r should equal the 
gravitational pull GM(r)/r 2 , and the radial dependence of v would then be expected to 
follow Kepler's law 

v 2 = GM(r)/r. (12) 

The surprising result for spiral galaxy rotation curves is, that the velocity does not 
follow Kepler's inverse-root law, but stays rather constant after attaining a maximum. 
The most obvious solution to this is that the galaxies are embedded in extensive, diffuse 
halos of dark matter If the mass M(r) enclosed inside the radius r, is proportional to r 
it follows that v(r) ~ constant. 

The rotation curve of most galaxies can be fitted by the superposition of 
contributions from the stellar and gaseous disks, sometimes a bulge, and the dark halo, 
modeled by a quasi - isothermal sphere. The inner part is difficult to model because 
the density of stars is high, rendering observations of individual star velocities difficult. 
Thus the fits are not unique, the relative contributions of disk and dark matter halo is 
model-dependent, and it is sometimes not even sure whether galactic disks do contain 
dark matter. Typically, dark matter constitutes about half of the total mass. 

In Fig. 3 we show the rotation curves fitted for eleven well-measured galaxies (25] 
of increasing halo mass. One notes, that the central dark halo component is indeed 
much smaller than the luminous disk component. At large radii, however, the need for 
a DM halo is obvious. On galactic scales, the contribution of DM generally dominates 
the total mass. Note the contribution of the baryonic component, negligible for light 
masses but increasingly important in the larger structures. 

The mass discrepancy emerges also as a disagreement between light and mass 
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Figure 3. Best disk - halo fits to the Universal Rotation Curve (dotted line is 
disk, dashed line is halo). Each object is identified by the halo virial mass, increasing 
downwards. From P. Salucci & al. [25] . 



distributions: light does not trace mass, the ratio 

(dM/dr)/(dL/dr) (13) 

is not constant, but increases with radius [26J. 

Gentile & al. [27] have shown that cusped profiles are in clear conflict with data on 
spiral galaxies. Central densities are rather flat, scaling approximately as p ^ luminous- 
The best-fit disk + NFW halo mass model fits the rotation curves poorly, it implies an 



Astrophysical and cosmological probes of dark matter 



11 




Figure 4. Left: The 4-dimensional space of luminosity, core radius, halo central 
density and fraction of DM. Right: The smooth surface of spiral galaxy rotation curves 
in the space of normalized radius R/R op ticai, magnitude M and rotation velocity V in 
km s _1 . P. Salucci priv. comm. and ref. |25j . 

implausibly low stellar mass-to-light ratio and an unphysically high halo mass. Clearly 
the actual profiles are of very uncertain origin. 

One notes in Fig. 3 that the shape of the rotation curve depends on the halo 
virial mass so that the distribution of gravitating matter, unlike luminous matter, is 
luminosity dependent. The old idea that the rotation curve stays constant after attaining 
a maximum is thus a simplification of the real situation. The rotation velocity can be 
expressed by a Universal Rotation Curve [25J. All spiral galaxies lie on a curve in the 
4-dimensional space of luminosity, core radius, halo central density and fraction of DM, 
see Fig. 4 

Our Galaxy is complicated because of what appears to be a noticeable density dip 
at 9 kpc and a smaller dip at 3 kpc, as is seen in Fig. 5 [28]. To fit the measured rotation 
curve one needs at least three contributing components: a central bulge, the star disk 
+ gas, and a DM halo [281 EU [30]. For small radii there is a choice of empirical rotation 
curves, and no DM component appears to be needed until radii beyond 15 kpc. 
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5. Strong and weak lensing 

A consequence of the Strong Equivalence Principle (SEP) is that a photon in a 
gravitational field moves as if it possessed mass, and light rays therefore bend around 
gravitating masses. Thus celestial bodies can serve as gravitational lenses probing the 
gravitational field, whether baryonic or dark without distinction. 

Since photons are neither emitted nor absorbed in the process of gravitational light 
deflection, the surface brightness of lensed sources remains unchanged. Changing the 
size of the cross-section of a light bundle only changes the flux observed from a source 
and magnifies it at fixed surface-brightness level. If the mass of the lensing object is very 
small, one will merely observe a magnification of the brightness of the lensed object an 
effect called microlensing. Microlensing of distant quasars by compact lensing objects 
(stars, planets) has also been observed and used for estimating the mass distribution of 
the lens-quasar systems. 

In Strong Lensing the photons move along geodesies in a strong gravitational 
potential which distorts space as well as time, causing larger deflection angles and 
requiring the full theory of General Relativity. The images in the observer plane can then 
become quite complicated because there may be more than one null geodesic connecting 
source and observer; it may not even be possible to find a unique mapping onto the source 
plane c/Fig.6 . Strong lensing is a tool for testing the distribution of mass in the lens 
rather than purely a tool for testing General Relativity. An illustration is seen in Fig. 
7 where the lens is an elliptical galaxy [32J . 
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Figure 6. Wave fronts and light rays in the presence of a cluster perturbation. From 
N. Straumann 31 . 



At cosmological distances one may observe lensing by composed objects such as 
galaxy groups which are ensembles of "point-like" , individual galaxies. Lensing effects 
are very model-dependent, so to learn the true magnification effect one needs very 
detailed information on the structure of the lens. 

Weak Lensing refers to deflection through a small angle when the light ray can be 
treated as a straight line (Fig. 6), and the deflection as if it occurred discontinuously 
at the point of closest approach (the thin-lens approximation in optics). One then only 
invokes SEP to account for the distortion of clock rates. 

The large-scale distribution of matter in the Universe is inhomogeneous in every 
direction, so one can expect that everything we observe is displaced and distorted by 
weak lensing. Since the tidal gravitational field and the deflection angles depend neither 
on the nature of the matter nor on its physical state, light deflection probes the total 
projected mass distribution. Lensing in infrared light offers an additional advantage of 
being able to sense distant background galaxies, since their number density is higher 
than in the optical range. 

Background galaxies would be ideal tracers of distortions if they were intrinsically 
circular, because lensing transforms circular sources into ellipses. Any measured 
ellipticity would then directly reflect the action of the gravitational tidal field of the 
interposed lensing matter, and the statistical properties of the distortions would reflect 
the properties of the matter distribution. But many galaxies are actually intrinsically 
elliptical, and the ellipses are randomly oriented. This introduces noise into the inference 
of the tidal field from observed ellipticities. A useful feature in the sky is a fine-grained 
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Figure 7. This image resulted from color-subtraction of a lensing singular isothermal 
elliptical galaxy. The strongly lensed object forms two prominent arcs A, B and a less 
extended third image C. From R.J. Smith & al. [32] 

pattern of faint and distant blue galaxies appearing as a 'wall paper'. This makes 
statistical weak-lensing studies possible, because it allows the detection of the coherent 
distortions imprinted by gravitational lensing on the images of the galaxy population. 

Thus weak lensing has become an important technique to map non-luminous 
matter. A reconstruction of one of the largest and most detailed weak lensing surveys 
undertaken with the Hubble Space Telescope is shown in Fig. 8 [33]. This map covers 
a large enough area to see extended filamentary structures. 

A very large review on lensing by R. Massey & al. [21] can be recommended. We 
show several examples of lensing by clusters in Sec. 15 

6. Elliptical galaxies 

Elliptical galaxies are quite compact objects which mostly do not rotate so their mass 
cannot be derived from rotation curves. The total dynamical mass is then the virial 
mass as derived from the velocity dispersions of stars and the anisotropies of their orbits. 
However, to disentangle the total mass profile into its dark and its stellar components 
is not straightforward, because the dynamical mass decomposition of dispersions is 
not unique. The luminous matter in the form of visible stars is a crucial quantity, 
indispensable to infer the dark component. When available one also makes use of strong 
and weak lensing data, and of the X-ray properties of the emitting hot gas. The gravity 
is then balanced by pressure gradients as given by Jeans' Equation. 

Inside the half light radius R e the contribution of the dark matter halo to the central 
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Figure 8. Map of the dark matter distribution in the 2-square degree COSMOS field: 
the linear blue scale on top shows the gravitational lensing magnification k, which is 
proportional to the projected mass along the line of sight. From R. Massey & al. [33J 

velocity dispersion is often very small, < 100 km s -1 , so that the dark matter profile is 
intrinsically unresolvable. The outer mass profile is compatible with NFW, Eq.(j7l), and 
with Burkert, Eq.Q, as well. Important information on the mass distribution can be 
obtained from the Fundamental Plane, Eq. (|6]). which yields the coefficients a = 1.8, 
b = 0.8. Note that this is in some tension with the Virial Theorem, perhaps due to 
variations in the central dispersions, cr , of the stellar populations. 

O. Tiret & al. [35] concluded from a study of 23 giant elliptical galaxies with central 
velocity dispersions > 330 km s -1 , that the mass within 5 — 10 kpc is dominated by the 
stars, not by DM. On the average the dark matter component contributes less than 5% 
to the total velocity dispersions. 

The ELIXR survey is a volume-limited (< 110 Mpc) study by P. J. Humphrey 
& al. [55] . of optically selected, isolated, L* elliptical galaxies in particular the NGC 
1521, for which X-ray data from Chandra and XMM exist. The isolation condition 
selects the appropriate galaxy halo and reduces the influence of a possible group-scale 
or cluster-scale halo. 

Most of the baryons are in a morphologically relaxed hot gas halo detectable out 
to ~ 200 kpc, that is well described by hydrostatic models. The baryons and the dark 
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Figure 9. Radial mass profile of the elliptical galaxy NGC 1521 from a model 
calculation (not fitted to the measured points shown). The solid black line indicates 
the total enclosed mass (la errors in grey), the dashed red line is the stellar mass, the 
dotted blue line is the dark matter, and the dash-dot magenta line is the gas mass 
contribution. From P. J. Humphrey & al. [55] . 

matter conspire to produce a total mass density profile that can be well-approximated 
by a power law, p to t oc r~ a over a wide range (as has been noted before, see references 
in [351 [36]). 

The fitting method involves solving the equation of hydrostatic equilibrium to 
compute temperature and density profile models, given parametrized mass and entropy 
profiles. The models are then projected onto the sky and fitted to the projected 
temperature and density profiles. A fit ignoring DM was poor, but inclusion of DM 
improved the fit highly significantly: DM was required at 8.2cr. We show this fit in Fig. 
9. In several studies [351 f° r most of the radii the dark matter contribution is very 
small although statistically significant. 

7. Mass to luminosity ratios and dwarf spheroidals 

The mass-to-light ratio of an astronomical object is defined as T e M/L. Stellar 
populations exhibit values T = 1 — 10 in solar units, in the solar neighborhood 
T = 2.5 - 7, in the Galactic disk T = 1.0 - 1.7 from C. Flynn & al. [3E]. 
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Dwarf spheroidal galaxies (dSph) are the smallest stellar systems containing dark 
matter and exhibit very high M/L ratios, T = 10 — 100. In Andromeda IX T= 93 
+120/-50, in Draco T = 330 ± 125. The dwarf spheroidals have radii of ~ 100 pc 
and central velocity dispersions ~ 10 km s _1 which is larger than expected for self- 
gravitating, equilibrium stellar populations. The generally accepted picture has been, 
that dwarf galaxies have slowly rising rotation curves and are dominated by dark matter 
at all radii. 

However, R. A. Swaters & al. |39j have reported observations of H I rotation 
curves for a sample of 73 dwarf galaxies, among which eight galaxies have sufficiently 
extended rotation curves to permit reliable determination of the core radius and the 
central density. They found that dark matter only becomes important at radii larger 
than three or four disk scale lengths. Their conclusion is, that the stellar disk can 
explain the mass distribution over the optical parts of the galaxy, and dark matter only 
becomes relevant at large radii. However, the required stellar mass-to-light ratios are 
high, up to 15 in the R-band. 

Comparing the properties of dwarf galaxies in both the core and outskirts of the 
Perseus Cluster, Penny and Conselice (10] found a clear correlation between mass-to- 
light ratio and the luminosity of the dwarfs, such that the faintest dwarfs require the 
largest fractions of dark matter to remain bound. This is to be expected, as the fainter 
a galaxy is, the less luminous mass it will contain, therefore the higher its dark matter 
content must be to prevent its disruption. Dwarfs are more easily influenced by their 
environment than more massive galaxies 

The distance to the Perseus Cluster prevents an easy determination of T, so S. J. 
Penny & C. J. Conselice (10] instead determined the dark matter content of the dwarfs 
by calculating the minimum mass needed in order to prevent tidal disruption by the 
cluster potential, using their sizes, the projected distance from the cluster center to 
each dwarf and the mass of the cluster interior. Three of 15 dwarfs turned out to have 
mass-to-light ratios smaller than 3, indicating that they do not require dark matter. 

Ultra-compact dwarf galaxies (UCDs) are stellar systems with masses of around 
10 7 — 10 8 M sun and half-mass radii of 10-100 pc. A remarkable properties of UCDs is 
that their dynamical mass-to-light ratios are on average about twice as large as those 
of globular clusters of comparable metallicity, and also tend to be larger than what one 
would expect based on simple stellar evolution models. UCDs appear to contain very 
little or no dark matter. 

H. Baumgardt & S. Mieske jl2] have presented collisional N-body simulations which 
study the coevolution of a system composed of stars and dark matter. They find that 
DM gets removed from the central regions of such systems due to dynamical friction and 
mass segregation of stars. The friction timescale is significantly shorter than a Hubble 
time for typical globular clusters, while most UCDs have friction times much longer 
than a Hubble time. Therefore, a significant dark matter fraction may remain within 
the half-mass radius of present-day UCDs, making dark matter a viable explanation for 
their elevated mass-to-light ratios. 
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A different type of systems are the ultra-faint dwarf galaxies (UFDs). When 
interpreted as steady state objects in virial equilibrium by V. Belokurov & al. [41], 
they would be the most DM dominated objects known in the Universe. Their half-light 
radii range from 70 pc to 320 pc. 

A special case is the UFD disk galaxy Segue 1, studied by M. Xiang-Gruess & al. 
|43j . which has a baryon mass of only about 1000 solar masses. One interpretation is 
that this is a thin non-rotating stellar disk not accompanied by a gas disk, embedded 
in an axisymmetric DM halo and with a ratio / = Mh a io/Mb ~ 200. But if the disk 
rotates, / could be as high as 2000. If Segue 1 also has a magnetized gas disk, the dark 
matter halo has to confine the effective pressure in the stellar disk and the magnetic 
Lorentz force in the gas disk as well as possible rotation. Then / could be very large 
[43] . Another interpretation is that Segue 1 is an extended globular cluster rather than 
an UFD [41]. 

8. Small galaxy groups emitting X-rays 

There are examples of groups formed by a small number of galaxies which are enveloped 
in a large cloud of hot gas (ICM), visible by its X-ray emission. One may assume that 
the electron density distribution associated with the X-ray brightness is in hydrostatic 
equilibrium, and one can extract the ICM radial density profiles by fits. 

The amount of matter in the form of hot gas can be deduced from the intensity of 
this radiation. Adding the gas mass to the observed luminous matter, the total amount 
of baryonic matter, M&, can be estimated, see M. Markevitch & al. [HJ and C. De Boni 
& G. Bertin |45]. In clusters studied, the gas fraction increases with the distance from 
the center; the dark matter appears more concentrated than the visible matter. 

The temperature of the gas depends on the strength of the gravitational field, from 
which the total amount of gravitating matter, M grav , in the system can be deduced. In 
many such small galaxy groups one finds M gI3 , v /Mb > 3, testifying to a dark halo present. 
An accurate estimate of M grav requires that also dark energy is taken into account, 
because it reduces the strength of the gravitational potential. There are sometimes 
doubts whether all galaxies appearing near these groups are physical members. If not, 
they will artificially increase the velocity scatter and thus lead to larger virial masses. 

On the scale of large clusters of galaxies like the Coma, it is generally observed that 
DM represents about 85% of the total mass and that the visible matter is mostly in the 
form of a hot ICM. 

9. Mass autocorrelation functions 

If galaxy formation is a local process, then on large scales galaxies must trace mass. This 
requires the study of how galaxies populate DM halos. In simulations one attempts to 
track galaxy and DM halo evolution across cosmic time in a physically consistent way, 
providing positions, velocities, star formation histories and other physical properties for 
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Figure 10. Dark matter halo mass Mhaio as a function of stellar mass M*. The thick 
black curve is the prediction from abundance matching assuming no dispersion in the 
relation between the two masses. Red and green dashed curves assume some dispersion 
in logM*. The dashed black curve is the satellite fraction as a function of stellar mass, 
as labeled on the axis at the right-hand side of the plot. From Qi Guo & al. [15] , 



the galaxy populations of interest. 

Guo & al. [IB] use abundance matching arguments to derive an accurate relation 
between galaxy stellar mass and DM halo mass. They combine a stellar mass function 
based on spectroscopic observations with a precise halo/subhalo mass function obtained 
from simulations. Assuming this stellar mass - halo mass relation to be unique and 
monotonic, they compare it with direct observational estimates of the mean mass of 
halos surrounding galaxies of given stellar mass inferred from gravitational lensing and 
satellite galaxy dynamics data, and use it to populate halos in simulations. The stellar 
mass - halo mass relation is shown in Fig. 10. 

The implied spatial clustering of stellar mass turns out to be in remarkably 
good agreement with a direct and precise measurement. By comparing the galaxy 
autocorrelation function with the total mass autocorrelation function, as averaged over 
the Local Supercluster (LSC) volume, one concludes that a large amount of matter in 
the LSC is dark. 

A similar study is that of Boyarsky Sz al. (17] who find a universal relation between 
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Figure 11. Dark matter column density vs. dark matter halo mass in solar units. 
From A. Boyarsky & al. [47] 



DM column density and DM halo mass, satisfied by matter distributions at all observable 
scales in halo sizes from 10 8 to 10 16 M sun , as shown in Fig. 11. Such a universal property 
is difficult to explain without dark matter. 

10. Cosmic Microwave Background (CMB) 

The tight coupling between radiation and matter density before decoupling caused the 
primordial adiabatic perturbations to oscillate in phase. Beginning from the time of last 
scattering, the receding horizon has been revealing these frozen density perturbations, 
setting up a pattern of standing acoustic waves in the baryon-photon fluid. After 
decoupling, this pattern is visible today as temperature anisotropics with a certain 
regularity across the sky. 

The primordial photons are polarized by the anisotropic Thomson scattering 
process, but as long as the photons continue to meet free electrons their polarization is 
washed out, and no net polarization is produced. At a photon's last scattering however, 
the induced polarization remains and the subsequently free-streaming photon possesses 
a quadrupole moment. 

Temperature and polarization fluctuations are analyzed in terms of multipole 
components or powers. The resulting distribution of powers versus multipole £, or 
multipole moment k = 2ir/£, is the power spectrum which exhibits conspicuous Doppler 
peaks. In Fig. 12 we display the radiation temperature (TT) and temperature- re- 
polarization correlation (TE) power spectra from the 7-year data of WMAP as functions 
of multipole moments jlH] . The spectra can then be compared to theory, and theoretical 
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Figure 12. The CMB radiation temperature (TT) and temperature-polarization (TE) 
power spectra from the seven-year WMAP 94 GHz maps. show anisotropics which can 
be analyzed by power spectra as functions of multipole moments. The solid line shows 
the best-fit prediction for the flat ACDM model. From D. Larson & al. [4T?] 

parameters determined. Many experiments have determined the power spectra, so a 
wealth of data exists. 

Baryonic matter feels attractive self-gravity and is pressure-supported, whereas dark 
matter only feels attractive self-gravity, but is pressureless. Thus the Doppler peaks 
in the CMBR power spectrum testify about baryonic and DM, whereas the troughs 
testify about rarefaction caused by the baryonic pressure. The position of the first peak 
determines VL m h 2 . Combining the TT data with determinations of the Hubble constant 
h, the WMAP team can determine the total mass density parameter Q m = f2& + fldm- 
The ratio of amplitudes of the second-to-first Doppler peaks determines the baryonic 
density parameter to be fi& = 0.0449 ± 0.0028 and the dark matter component to be 
Q dm = 0.222 ± 0.026 gg, thus Q m = 0.267. 

Power spectra at higher multipole moments have been measured with the Atacama 
Cosmology Telescope (ACT) [50] at 148 GHz and 218 GHz, as well as the cross-frequency 
spectrum between these two channels, and found to be in agreement with the 7-year 
WMAP 94 GHz maps in the common range 400 < £ < 1000. The ACT has also been 
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Figure 13. The ACT 148GHz power spectrum multiplied by £ 4 is shown for lensed 
(orange curve) and unlensed models (green curve). From S. Das & al. 50J. 

able to measure the lensing of the CMB signal at a significance of 2.8a, which slightly 
smooths out the acoustic peaks, cf Fig. 13. 

In a fit of the flat ACDM model to the data the dark matter density parameter 
comes out slightly higher than WMAP and the baryonic density slightly lower so the 
total density parameter for WMAP and ACT added is Q m = 0.276 ± 0.016 [SI]. 
Information on the TE correlations comes from several measurements, among them 
WMAP [49], and on the E-mode polarization power spectrum alone (EE) from the 
QUAD collaboration [52], Fig. 14. 

The results show two surprises: Firstly, since fl m <C 1, a large component Q\ ~ 0.74 
is missing, of unknown nature, and termed dark energy. The second surprise is that 
ordinary baryonic matter is only a small fraction of the total matter budget. The 
remainder is then dark matter, of unknown composition. Of the 4.5% of baryons in the 
Universe only about 1% is stars. 

11. Baryonic acoustic oscillations (BAO) 

A cornerstone of cosmology is the Copernican principle, that matter in the Universe 
is distributed homogeneously, if only on the largest scales of superclusters separated 
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Figure 14. The E-mode polarization power spectrum (EE) from the CMB 
observations of the QUaD collaboration & al. [55] 



by voids. On smaller scales we observe inhomogeneities in the forms of galaxies, 
galaxy groups, and clusters. The common approach to this situation is to turn 
to non-relativistic hydrodynamics and treat matter in the Universe as an adiabatic, 
viscous, non-static fluid, in which random fluctuations around the mean density appear, 
manifested by compressions in some regions and rarefactions in other. The origin of 
these density fluctuations was the tight coupling established before decoupling between 
radiation and charged matter density, causing them to oscillate in phase. An ordinary 
fluid is dominated by the material pressure, but in the fluid of our Universe three effects 
are competing: gravitational attraction, density dilution due to the Hubble flow, and 
radiation pressure felt by charged particles only. 

The inflationary fluctuations crossed the post-inflationary Hubble radius, to come 
back into vision with a wavelength corresponding to the size of the Hubble radius at 
that moment. At time t eq the overdensities began to amplify and grow into larger 
inhomogeneities. In overdense regions where the gravitational forces dominate, matter 
contracts locally and attracts surrounding matter, becoming increasingly unstable until 
it eventually collapses into a gravitationally bound object. In regions where the pressure 
forces dominate, the fluctuations move with constant amplitude as sound waves in the 
fluid, transporting energy from one region of space to another. 

Inflationary models predict that the primordial mass density fluctuations should 
be adiabatic, Gaussian, and exhibit the same scale invariance as the CMB fluctuations. 
The baryonic acoustic oscillations can be treated similarly to CMB, they are specified 
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Figure 15. BAO in power spectra calculated from (a) the combined SDSS and 
2dFGRS main galaxies, (b) the SDSS DR5 LRG sample, and (c) the combination 
of these two samples (solid symbols with la errors). The data are correlated and the 
errors are calculated from the diagonal terms in the covariance matrix. A Standard 
ACDM distance-redshift relation was assumed to calculate the power spectra with 
fl m = 0.25, fl A = 0.75. From W. J. Percival & al. [53]. 



by the dimensionless mass autocorrelation function which is the Fourier transform of 
the power spectrum of a spherical harmonic expansion. The power spectrum is shown 
in Fig. 15 [53] . 

As the Universe approached decoupling, the photon mean free path increased 
and radiation could diffuse from overdense regions into underdense ones, thereby 
smoothing out any inhomogeneities in the plasma. The situation changed dramatically 
at recombination, at time 380 000 yr after Big Bang, when all the free electrons 
suddenly disappeared, captured into atomic Bohr orbits, and the radiation pressure 
almost vanished. Now the baryon acoustic waves and the CMB continued to oscillate 
independently, but adiabatically, and the density perturbations which had entered the 
Hubble radius since then could grow with full vigor into baryonic structures. 

The scale of BAO depends on Q m and on the Hubble constant, h, so one needs 
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information on h to break the degeneracy. The result is then Vt m f» 0.26. In the ratio 
Qt,/Q m the /i-dependence cancels out, so one can also quantify the amount of DM on 
very large scales by Qf,/Q m = 0.18 ± 0.04. 

12. Galaxy formation in purely baryonic matter? 

We have seen in Sec. 10 that the baryonic density parameter, f^> is very small. The 
critical density £l cr u is determined by the expansion speed of the Universe, and the 
mean baryonic density of the Universe (stars, interstellar and intergalactic gas) is only 

n b = o.o45 ng. 

The question arises whether the galaxies could have formed from primordial density 
fluctuations in a purely baryonic medium. We have also noted, that the fluctuations in 
CMB and BAO maintain adiabaticity. The amplitude of the primordial baryon density 
fluctuations would have needed to be very large in order to form the observed number 
of galaxies. But then the amplitude of the CMB fluctuations would also have been very 
large, leading to intolerably large CMB anisotropies today. Thus galaxy formation in 
purely baryonic matter is ruled out by this argument alone. 

Thus one concludes, that the galaxies could only have been formed in the presence 
of gravitating dark matter which started to fluctuate early, unhindered by radiation 
pressure. This conclusion is further strengthened in the next Section. 

13. Large Scale Structures simulated 

In the ACDM paradigm, the nonlinear growth of DM structure is a well-posed problem 
where both the initial conditions and the evolution equations are known (at least when 
the effects of the baryons can be neglected). 

The Aquarius Project |54j is a Virgo Consortium program to carry out high- 
resolution DM simulations of Milky- Way-sized halos in the ACDM cosmology. This 
project seeks clues to the formation of galaxies and to the nature of the dark matter by 
designing strategies for exploring the formation of our Galaxy and its luminous and dark 
satellites. The galaxy population on scales from 50 kpc to the size of the observable 
Universe has been predicted by hierarchical ACDM scenarios, and compared directly 
with a wide array of observations. So far, the ACDM paradigm has passed these tests 
successfully, particularly those that consider the large-scale matter distribution and has 
led to the discovery of a universal internal structure for DM halos. As was noted in Sec. 
12, the observed structure of galaxies, clusters and superclusters, as illustrated by Fig. 
16, could not have formed in a baryonic medium devoid of dark matter. 

Given this success, it is important to test ACDM predictions also on smaller scales, 
not least because these are sensitive to the nature of the dark matter. Indeed, a 
number of serious challenges to the paradigm have emerged on the scale of individual 
galaxies and their central structure. In particular, the abundance of small DM subhalos 
predicted within CDM halos is much larger than the number of known satellite galaxies 
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Figure 16. The left panel shows the projected dark matter density at z = in a 
slice of thickness 13.7 Mpc through the full box (137 Mpc on a side) of the 900 3 
parent simulation. The right panel show this halo resimulated at a different numerical 
resolution. The image brightness is proportional to the logarithm of the squared DM 
density projected along the line-of-sight. The circles mark r 50 , the radius within which 
the mean density is 50 times the background density. From V. Springel & al. [53] 

surrounding the Milky Way (M. Boylan-Konchin & al. [IH] and references therein). 
14. Dark matter from overall fits 

In Sec. 10 we have seen that the WMAP 7-year CMB data together with the Hubble 
constant value testify about the existence of DM [39], [51]. In Sec. 11 we addressed 
the BAO data [53J with the same conclusion. In overall fits one combines these with 
supernova data (SN la) which offer a constraint nearly orthogonal to that of CMB in the 
— fi m -plane. The Union compilation of 307 selected SN la includes the recent large 
samples of SNe la from the Supernova Legacy Survey, the ESSENCE Survey, the older 
data sets, as well as the recently extended data set of distant supernovae observed with 
HST. M. Kowalski & al. [55] present the latest results from this compilation and discuss 
the cosmological constraints and its combination with CMB and BAO measurements. 
The CMB constraint is close to the line f^A + — 1> whereas the supernova constraint 
is close to the line Q\ — 1.6 x Q m = 0.2. The BAO data constrain Q m , but hardly at all 
f^A- This is shown in Fig. 17. 

Defining the vacuum energy density parameter by Qk = 1 — a flat Universe 

corresponds to Qk = 0. For a non-flat AC DM Universe with a cosmological constant 
responsible for dark energy, a simultaneous fit to the data sets gives 

tt m = 0.285+0.020/-0.019±0.011,n fc = -0.009+0.009+0.002/-0.010-0.003 

where the first error is statistical and the second error systematic. Clearly one notes 
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Figure 17. 68.3 %, 95.4 % and 99.7% confidence level contours on and fl m obtained 
from CMB, BAO and the Union SN set, as well as their combination (assuming w = 
-1). Note the straight line corresponding to a flat Universe with fl\ + Q m — 1. From 
M. Kowalski & al. [S5]. 
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that the Universe is consistent with being flat. Subtracting fi& = 0.045 from Vt m = 0.285 
one obtains the density parameter for DM, Q^m ~ 0.24. Assuming flatness, M. Kowalski 
& al. (55] find Q m = 0.274 ±0.016 ±0.013. This compares well with the combined 7-year 
WMAP data and the ACT data, fi m = 0.276 ± 0.016 [SI]. If one fits different models 
having more free parameters, one gets slightly different results, but all within these la 
errors. 
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Figure 18. The merging cluster 1E0657-558. On the right is the smaller Bullet cluster 
which has traversed the larger cluster. The colors indicate the X-ray temperature of 
the plasma: blue is coolest and white is hottest. The green contours are the weak 
lensing reconstruction of the gravitational potential of the cluster. From D. Clowe & 
al. [55] 

15. Merging galaxy clusters 

In isolated galaxies and galaxy clusters all matter components contributing to the 
common gravitational potential are more or less centrally-symmetrically coincident. 
This makes the dissection of DM from the baryonic components difficult and dependent 
on parametrization, as we have discussed in Sec. 3. In merging galaxy clusters however, 
the separate distributions of galaxies, intracluster gas and DM may become spatially 
segregated permitting separate observations. The visually observable galaxies behave 
as collisionless particles, the baryonic intracluster plasma is fluid-like, experiences ram 
pressure and emits X-rays, but non-interacting DM does not feel that pressure, it only 
makes itself felt by its contribution to the common gravitational potential. 

Major cluster mergers are the most energetic events in the Universe since the Big 
Bang. Shock fronts in the intracluster gas are the key observational tools in the study 
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of these systems. When a subcluster traverses a larger cluster it cannot be treated as a 
solid body with constant mass moving at constant velocity. During its passage through 
the gravitational potential of the main cluster it is shrinking over time, stripped of gas 
envelope and decelerating. Depending on the ratio of the cluster masses, the gas forms 
a bow shock in front of the main cluster, and this can even be reversed at the time when 
the potentials coincide. 

We shall now meet several examples of galaxy cluster mergers where the presence 
of DM could be inferred from the separation of the gravitational potential from the 
position of the radiating plasma. 

15.1. The Bullet cluster 1E0657-558 

The exceptionally hot and X-ray luminous galaxy cluster 1E0657-558, the Bullet cluster 
at redshift z = 0.296, was discovered by Tucker et al. in 1995 [56j in Chandra X-ray 
data. Its structure as a merger of a 2.3 xlO 14 M sun subcluster with a main 2.8 xlO 14 
M sun cluster was demonstrated by Markevitch et al. EE] and Clowe et al. [60, 59J. 
This was presented as the first clear example of a bow shock in a heated intracluster 
plasma. 

With the advent of high-resolution lensing Bradac et al. [611 E2] developed a 
technique combining multiple strongly- lensed Hubble Space Telescope multi-color images 
of identified galaxies, with weakly lensed and elliptically distorted background sources. 
The reconstructed gravitational potential does not trace the X-ray plasma distribution 
which is the dominant baryonic mass component, but rather approximately traces the 
distribution of bright cluster member galaxies, cf Fig. 18. 

The center of the total mass is offset from the center of the baryonic mass peaks, 
proving that the majority of the matter in the system is unseen. In front of the bullet 
cluster which has traversed the larger one about 100 Myr ago with a relative velocity 
of 4500 km s -1 , a bow shock is evident in the X-rays. The main cluster peak and the 
distinct subcluster mass concentration are both clearly offset from the location of the 
X-ray gas [62] . A recent analysis of this system [72J confirms the results of references 
\59\ ISTH 162] , and in addition finds that dark matter forms three distinct clumps. 

15.2. The galaxy cluster pair MACS J0025. 4-1222 

Another merging system with similar characteristics but with lower spatial resolution has 
been reported by Bradac et al. [19J, the post-merging galaxy cluster pair MACS J0025.4- 
1222, also called the Baby Bullet. It has an apparently simple geometry, consisting of two 
large subclusters of similar richness, about 2.5 xlO 14 M sun , both at redshift z = 0.586, 
colliding in approximately the plane of the sky. Multiple images due to strong lensing 
of four distinct components could be identified. The combined strong and weak lensing 
analysis follows the method in ref. [62J. 




The two distinct mass peaks are clearly offset by 4cr from the main baryonic 
component, which is the radiating hot gas observed by Chandra. The relative merging 




Astrophysical and cosmological probes of dark matter 



30 




Figure 19. The color composite of the cluster MACS J0025.4-1222. Overlaid in 
red contours is the surface mass density (linearly spaced) from the combined weak 
and strong lensing mass reconstruction. The X-ray brightness contours (also linearly 
spaced) are overlaid in yellow and the I-band light is overlaid in white. The measured 
peak positions and error bars for the total mass of the two cluster components are 
shown as cyan crosses. From M. Bradac & al. [15] . 

velocity is estimated to be 2000 km s™ 1 . In Fig. 19 we show linearly spaced surface mass 
density contours and X-ray brightness contours. The majority of the mass is spatially 
coincident with the identified galaxies which implies, that the cluster must be dominated 
by a relatively collisionless form of dark matter. 

15.3. The merging system A 1758 

A much more complicated merging system is A1758 at redshift z = 0.279, analyzed 
by the same team as above, B. Ragozzine & al. [53], and consisting of four clusters 
undergoing two separate mergers. The weak lensing mass peaks of the two northern 
clusters A1758N are separated at the 2.5 a level, whereas the two southern clusters are 
not well separated and have a disturbed X-ray morphology. There is no evidence for a 
merger between A1758N and A1758S in the X-ray signature and they have a projected 
separation of 2.0 Mpc. Note however the SZ results from the Arcminute Microkelvin 
Imager (AMI) in Cambridge (UK) on this system [Ml 165] . which sees a hint of a signal 
between the A1758N and A1758S. 
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Figure 20. The A1758N merger from B. Ragozzine & al. [63]. The blue contours 
represent the weak lensing mass reconstruction made from a background galaxy density 
of 24.0 galaxies/arcmin 2 . The outer blue contour begins at surface mass density 
k = 0.07 and each contour increases in steps of 0.045 up to n — 0.34. The red contours 
follow the X-ray gas mass obtained in the Chandra exposure. The NW cluster's BCG 
aligns with the X-ray gas and the weak lensing peak. The SE cluster's BCG and weak 
lensing peak are well separated from the X-ray gas, which has a bright peak near the 
midpoint of the two weak lensing peaks. 



A1758N introduces a new geometry that is different from the previously discussed 
mergers: one weak lensing peak overlaps an X-ray peak, while the other weak lensing 
peak is clearly separated from the X-ray component, cf Fig. 20. 

Since no strong lensing has yet been confirmed, conclusions about cluster masses 
and DM would have to wait for better lensing data. 

15. 4- The merging cluster Abell 2146 

Chandra observations of the cluster Abell 2146 at a redshift of z — 0.234 have revealed 
two shock fronts, H. R. Russell & al. [66]. The X-ray morphology suggests a recent 
merger where a subcluster containing a dense core has passed through the center of a 
second cluster, the remnant of which appears as the concentration of gas to the NW. The 
strongly peaked core has just emerged from the primary core, and is trailing material 
that has been ram pressure stripped in the gravitational potential. This material appears 
as a warmer stream of gas behind the subcluster core, and trails back to the hottest 
region of the disrupted main cluster. 
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Figure 21. Left: Projected emission measure per unit area map for the merging 
cluster Abell 2146. Center: Projected temperature map. Right: Projected "pressure" 
map produced by multiplying the emission measure and temperature maps. From H. 
R. Russell & oL (SB] 

Four steep surface brightness edges can be defined: two in the SE sector in front of 
the subcluster core and another two in the NW sector, cf. Fig. 21. The interpretation 
is [HE] that an upstream shock is generated as the gravitational potential minimum 
fluctuates rapidly during the core passage, reaching an extreme minimum when the 
two cluster cores coalesce. This causes a significant amount of the outer cluster gas to 
flow inwards. When the subcluster core exits the main core the gravitational potential 
rapidly returns to its premerger level and expels much of the newly arrived gas which in 
turn collides with the residual infall, forming an inward traveling shock front. Behind 
the subcluster, the ambient cluster gas that was pushed aside during its passage will 
fall back and produce tail shocks. 

Since no weak lensing analysis is available as yet, nothing can be said about the 
possible role of collisionless dark matter. 

15.5. The merging cluster Abell 2744 

Newly acquired data with the Advanced Camera for Surveys on the Hubble Space 
Telescope, HST, shows that the cluster Abell 2744 is a complicated merger between 
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Figure 22. The surface-mass density contours of the merging cluster Abell 2744 are 
shown in cyan and the X-ray luminosity contours in magenta. The peak positions of 
the Southern core, the N, NW, and W clumps are indicated by the green likelihood 
contours. The small red circles show the positions of the local overdensities in the gas 
distribution, associated with each DM clump. The white rulers show the separation 
between DM peaks and the bright cluster galaxies and local gas peaks. From J. Merten 
& al. [67]. 

three or four separate bodies, as analyzed by J. Merten & al. [67J. The position 
and mass distribution of the Southern core have been tightly constrained by the strong 
lensing of 11 background galaxies producing 31 multiple images. The N and NW clumps 
lack such images from strong lensing, indicating that they are less massive. There is 
also weak lensing information from HST, VLT, and Subaru available. 

The joint gravitational lensing analysis combines all the strongly lensed multiply- 
imaged systems and their redshifts with weak lensing shear catalogues from all three 
telescopes to reconstruct the cluster's lensing potential, shown in Fig. 22. The Core, 
NW and W clumps are clear detections in the surface-mass density distribution with 
Ilex, 4.9a and 3.8a significance over background, respectively. Somewhat fainter with 
2.3cr significance is the N structure, but it clearly coincides with a prominent X-ray 
substructure found by M. S. Owers & al. [68 J. 

To determine the geometric configuration of the collision, the location of shock 
fronts and velocities, densities and temperatures in the intracluster medium, all existing 
X-ray data from Chandra [HE] were included and reanalyzed. Overlaying the lensing 
mass reconstruction and the luminosity contours of the emission in Fig. 23 shows 
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an extremely complex picture of separations between the dark matter and baryonic 
components. 

In the core region which is the most massive structure within the merging system, 
there is no large separation between the distributions of total mass and baryons. The 
separation of the peaks in the lensing and X-ray maps is similar to that in the Bullet 
cluster [HE] and Baby Bullet [19] . The Northern mass substructure is ~ 2.6 times lighter 
than the Core, and the X-ray emission lags behind the dark matter to the South. 
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Figure 23. The proposed merger scenario of the cluster Abell 2744 in time-ordered 
sequence. The NE-SW subclusters merge first (1) with the core, followed very soon 
(2) by the second merger, in the SE-NW direction. The gas slingshots (3) away to its 
present position at the extreme NW. In (4) we see the present setup. From J. Merten 
& al. [57]. 

The substructure in the Northwest is the second most massive and there might 
also be a second peak in the more Western area of the NW mass clump. However, it 
is difficult to ascertain whether this is a single, separate DM structure and to derive 
decisive separation between DM, X-ray luminous gas and bright cluster member galaxies. 
The X-ray peak to the Northwest of the NW2 mass peak appears to be an X-ray feature 
with no associated matter or galaxies, a "ghost" clump. 

One possible interpretation [67] of the complex merging scenario that has taken 
place in Abell 2744 is a near simultaneous double merger 0.12 — 0.15 Gyr ago. first in 
the NE-SW direction, cf. Fig. 23. The Western clump probably passed closest through 
the main cluster, as it had its ICM ram-pressure stripped completely. The second 
merger, in the SE-NW direction, could even have consisted of two small clumps falling 
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into the core, attracted by the core and the Northern and Western clumps. After a 
first core passage, gas initially trails its associated DM but, while the dark matter slows 
down, the gas slingshots past it due to a combination of low ram-pressure stripping and 
adiabatic expansion and cooling [68] , ending up as the "ghost" clump. This scenario 
still requires further observations as well as verification via numerical simulations. 




Figure 24. Isopleth contours in yellow of the number density of galaxies in the 
cluster ATC-CL J0102-4915. The white contours show the X-ray emission. From F. 
Menenteau fc al. [70] 

15.6. "El Gordo", the fat cluster ACT-CL J0102-4915 

The Atacama Cosmology Telescope has presented properties for an exceptionally 
massive merging cluster, the ACT-CL J0102-4915 nicknamed El Gordo at redshift 
z = 0.87. It was discovered by Marriage & al. [69] selected by its bright Sunyaev- 
Zel'dovich (SZ) effect, confirmed optically and through its Chandra X-ray data [TO]- It 
is the most significant SZ cluster detection to date by nearly a factor of two, with an 
SZ decrement comparable to the Bullet cluster 1E0657-558 [62J. 

As can be seen from Fig. 24, the galaxy distribution is double peaked, whereas 
the peak in the X-ray emission lies between the density peaks. The X-ray peak forms a 
relatively cool bullet of low entropy gas like in the 1E0657-558. The steep fall-off in the 
X-ray surface brightness towards the SE, as well as the "wake" in the main cluster gas 
toward the NW, indicate that the bullet is apparently moving toward the SE. The SZ 
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and X-ray peaks are offset similar to that reported for the bullet-like cluster Abell 2146 

In the absence of a weak lensing mass reconstruction, the galaxy distribution can 
only be used as a proxy for the total mass distribution. Thus to conclude that an offset 
between baryonic and DM has been demonstrated is yet premature. 

15.7. The cluster merger DLSCL J0916. 2+2951 

A newly discovered [7T] major cluster merger at z = 0.53 is DLSCL J0916. 2+2951, in 
which the collisional cluster gas has become clearly dissociated from the collisionless 
galaxies and dark matter. The cluster was identified using optical and weak-lensing 
observations as part of the Deep Lens Survey. Follow-up observations with Keck, Subaru, 
Hubble Space Telescope, and Chandra show that the cluster is a dissociative merger 
which constrain the DM self-interaction cross-section to am^ M < 7 cm 2 /g. The system 
is observed at least 0.7 ± 0.2 Gyr since first pass-through, thus providing a picture of 
cluster mergers 2-5 times further progressed than similar systems observed to date. 

16. Comments and conclusions 

What we have termed "dark matter" is generic for observed gravitational effects on 
all scales: galaxies, small and large galaxy groups, clusters and superclusters, CMB 
anisotropies over the full horizon, baryonic oscillations over large scales, and cosmic 
shear in the large-scale matter distribution. The correct explanation or nature of dark 
matter is not known, whether it implies unconventional particles or modifications to 
gravitational theory, but gravitational effects have convincingly proved its existence in 
some form. 

The few per cent of the mass of the Universe found as baryonic matter in stars and 
dust clouds is well accounted for by nucleosynthesis. If there exist particles which were 
very slow at time t eq when galaxy formation started, they could be candidates for cold 
dark matter. They must have become non-relativistic much earlier than the leptons, 
and then decoupled from the hot plasma. 

Whenever laboratory searches discover a new particle, it must pass several tests 
in order to be considered a viable DM candidate: it must be neutral, compatible 
with constraints on self-interactions (essentially collisionless), consistent with Big Bang 
nucleosynthesis, and match the appropriate relic density. It must be consistent 
with direct DM searches and gamma-ray constraints, it must leave stellar evolution 
unchanged, and be compatible with other astrophysical bounds. 

The total dynamical mass of an astronomical system is derivable from the velocity 
dispersions or the rotation velocities of its components via the use of the Virial Theorem 
or Kepler's law, respectively. A most important probe is strong gravitational lensing 
which measures the total mass, but also weak lensing, the oscillations in the Cosmic 
Microwave Background and in the ambient baryonic medium. Probes separating dark 
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matter from total matter require in addition observations of visible light, infrared 
radiation, X-rays, the Sunyaev-Zel'dovich effect, and supernovae. Depending on the 
system under study there are many ways to combine these tools using empirical halo 
models, simulating stellar population models and galaxy formation models, comparing 
mass-to-light ratios and mass autocorrelation functions. The most remarkable systems 
are merging galaxy clusters which, by their motion, separate non-collisional dark matter 
from optically visible galaxies and hot, radiating gas. 

Regardless of the nature of dark matter, all theories attempting to explain it share 
the burden to explain the gravitational effects described in here. Thus there remains 
much to be done. 
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